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Abstract. The multimode conductance of a closed ring is found within the framework of 
a scattering approach. The expression can be regarded as a generalization of the Landauer 
formula. The treatment is essentially classical because we assume short coherence time. Our 
starting point is the Kubo formalism, but we also use a master equation approach for the 
derivation. As an example we calculate the conductance of a multimode waveguide with an 
attached cavity. 



1. Introduction 

The notion of conductance has gone several transformation in the last century. In 
the mesoscopic community ^ following Landauer, it is customary nowadays to 
consider the open geometry that is described in Fig. la, where a device is attached 
to left and right reservoirs, and the bias is understood as emerging from a chemical 
potential difference. For a single mode device it is argued that the conductance is 
essentially the transmission, while for a multimode device with "spinless" electrons 



where g^ m is the transmission matrix. We can optionally assume that the chemical 
potential of the two reservoirs is the same, and consider the effect of an electro motive 
force (EMF) such that the voltage drop is concentrated across a segment of the device. 
More generally we can consider the problem of driving a current by changing the 
electrical potential at some region of the device. The latter more general problem 
is known as "quantum pumping" . The calculation of the "geometric" conductance 
in the latter case leads to the Buttiker-Pretre-Thomas (BPT) formula [3], which is a 
generalization of the Landauer formula Eq. (JTJ . 

It is quite natural to ask what happens if the two leads are detached from the 
reservoirs, and the system is closed into a ring as in Fig. lb 01 El- Still we can induce 
EMF by changing in time an Aharonov-Bohm magnetic flux, or we can change the 
potential in some region of the device, so as to get an electrical current. In spite 
of much interest in closed mesoscopic rings a straightforward answer to this simple 
question has not been given. We shall review later the main published statements 



In this paper we are interested in circumstances such that the leading result for 
the conductance is of classical nature. This is completely analogous to the discussion 
of diffusive rings in circumstances such that the leading result is given by the Drudc 
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Figure 1. Panel (a) displays the standard Landauer (open) geometry of connecting a conductor 
to a left and a right reservoirs. In this illustration the conductor is a rectangular waveguide to 
which a cavity is attached. In panel (b) the leads are joined together so as to form a ring. The 
motion in the ring is assumed to be chaotic due to the scattering in and out of the cavity. Panel 
(c) is the schematic electrical engineering representation of the system. In panel (d) the system 
is modeled as a network. The scattering region is described by the transition matrix g^. In 
(c) and in (d) we indicate the presence of the environment by the gray shading. In fact (a) can 
be regarded as a special case of (b) provided one assumes that the effect of the environment 
is to randomize the velocity within the wire region. The current is measured via the section 
x = x\. The Electro motive force (EMF) is realized by time dependent Aharonov Bohm flux. 
The voltage drop can be concentrated anywhere along the ring (say across x = xq). Setting a 
chemical potential difference in the setup of panel (a) can be regarded as such particular option. 

formula. We are going to assume that the coherence time is shorter than the time that 
it takes for an electron to encircle the ring. Thus, as far as the dynamics is concerned, 
our analysis is essentially of classical nature. The case of a fully coherent ring requires 
further analysis which is beyond the scope of the present work. Therefore it will be 
discussed in separate publications: The one-mode case in an d the multi-mode case 
in Q2]. 

It is important to realize that our assumption of short coherence time is plays 
analogous role to that which is played by the reservoirs in the Landauer-BPT 
formalism. In fact we are going to explain that the Landauer-BPT formalism can 
be regarded as a special limit of the problem that we are going to consider. This will 
further illuminate the formal discussion in 0]. 
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1.1. Motivation 

The interest in the response of small mesoscopic rings is long standing [51 El El 03 EH • 
Measurements of the the conductance of closed mesoscopic rings has been performed 
already 10 years ago In a practical experiment a large array of two dimensional 
rings is fabricated. The conductance measurement can be achieved via coupling to 
a highly sensitive electromagnetic superconducting micro-resonator. In such setup 
the EMF is realized by creating a current through a "wire" that spirals on top of 
the array, and the conductance of the rings is determined via their influence on the 
electrical circuit. Another possibility is to extract the conductance from the rate of 
Joule heating. The later can be deduced from a temperature difference measurement 
assuming that the thermal conductance is known. 

We do not think that there is any problem to design rings of the type which is 
illustrated in Fig. 1. The question is why to bother? This brings us to the theoretical 
motivation for this work. Past theoretical studies are quite well summarized by the 
paper "(Almost) everything you always wanted to know about the conductance of 
mesoscopic systems" |l()j and see references therein. The major interest was in diffusive 
rings, and the main issue was weak localization corrections to the Drude result taking 
into account the level statistics, the type of occupation, etc. Our interest is in a 
different type of configuration, which is motivated by "quantum chaos" studies. Our 
claim (see next subsection) that the conductance of a closed ring can be larger than 
the number of open modes is quite unorthodox. 

More importantly: Our approach to the problem of mesoscopic conductance has a 
practical appeal. The reason for the popularity of the scattering approach in mesoscopic 
physics is its "plug and play" feature. The experimentalist is able to characterize the 
scattering properties of his/her device, and then he/she is able to make a prediction 
regarding the conductance. It is only natural to extend this " plug and play" approach 
to the analysis of conductance of closed rings. This extension is far from being trivial, 
as we explain later (in section 6). 

We care to make clear distinction between classical and quantal effects. This 
helps to develop a better intuition for the physics of such devices. The approach 
in the present work is in the spirit of the Boltzmann picture: The scattering "cross 
section", which is possibly of quantum mechanical nature, is taken as an input, while 
the overall dynamics is assumed to be of classical nature. The implications of quantum 
interference will be discussed in future works |14l ITS] . 



1.2. specific results 



The conductance of a single mode ring (A4 
transmission gx is given by the expression 
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We shall show that the multimode generalization of this formula is: 
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where g^ m and g^ m are the MxM transmission and reflection blocks of the transition 
matrix, As an example we analyze the system of Fig. la. We find that the Landauer 
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formula gives 
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where a is the dimensionless size of the opening to the cavity. In contrast to that for 
the multi-mode conductance of the corresponding ring structure (Fig. lb) we get 



G 



e 



id 

a 



\n{2M)) 



M 



(•5) 



Unlike the case of the Landauer conductance, the result does not reflect the number 
of open modes. This is because the contribution of the low modes is singular in the 
limit of small a. Furthermore, the conductivity (conductance per channel) diverges 
logarithmically in the classical limit. 



1.3. Outline 

We define the model system in sections 2 and 3 and the notion of conductance in 
section 4. In section 5 we argue that the notion of conductance is meaningful even in 
the absence of a bath. We further discuss the role of the environment in section 6, 
where we make distinguish between various type of bath induced effects. 

As explained in section 7 the purpose of the "linear response analysis" is to find 
the stationary-like state of the driven system. The procedure is to assume that in the 
absence of driving the system would be in a (strict) stationary state, which we regard 
as the zero order solution. Then we try to find a first order solution (in the EMF) to 
the time dependent problem. 

In sections 8 and 9 we discuss the Kubo approach to linear response. We take 
the simplest route following Refs. |12l 1161 1171 HH| leading to the fluctuation-dissipation 
version of the Kubo formula. The application of the Kubo formula to the analysis of 
the single mode conductance is presented in sections 10 and 11. The relation to the 
Landauer result is clarified in section 12. 

There are many (equivalent) ways to do "linear response analysis" . It turns out 
that the derivation of the multi-mode conductance formula becomes more transparent 
by adopting a master equation approach. This is carried out in section 13. 

Throughout the paper, and in particular in the concluding section, we emphasize 
that our approach has a straightforward extension to the analysis of quantum pumping. 
The Kubo approach allows a better understanding of the role which is played by the 
environment, and makes it possible to bridge between the strict quantum adiabatic 
limit and the other extreme limit of having an open geometry. 



2. Setting up the model 

Consider non- interacting spinless electrons in a ring, as in Fig. lb. The one-particle 
Hamiltonian is 

H= -(p-e$a(r)) 2 + V(r) (6) 
m 

where m and e are the mass and the charge respectively. The vector potential which 
is associated with the flux 4> is described by 

ja{r)-dr = l (7) 
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The dimensionality of the ring is d. The ring consists of a "wire" region and a scattering 
region. The motion of the particle inside the ring is assumed to be globally chaotic. 
The coordinate along the wire will be denoted as x. The scattering region is located 
at x <~ 0. 

In the geometry of Fig. lb the "wire" is a d = 2 waveguide of width L±. Later 
we describe the waveguide as a set of d — 1 wires (Fig. Id) such that each "wire" 
corresponds to an open mode of the waveguide. The length of the ring L is assumed 
to be large compared with the scattering region so as to allow meaningful definition 
of a scattering matrix in the quantum mechanical analysis (evanescent modes are 
ignored) . 

The ring is driven by a time dependent Aharonov-Bohm flux. The EMF — $ is 
assumed to be constant. There are various ways to introduce the EMF into the ring. 
One possibility is to have all the voltage drop over a section at x — x . Namely, 

a(x) = 5(x — x ) (8) 

For sake of later analysis we define a generalized force which is associated with the 
flux: 

dl~t 

:F =--Q^ = e v5(x-x ) (9) 

where v is the velocity in the x direction. In the quantum mechanical case a 
symmetrization is implicit. This is in fact a current operator. Obviously we do not 
have to measure the current at the same point where we apply the voltage. So for 
sake of generality we introduce the notation 

l = ev8{x-x l ) (10) 

We also note that with uniform averaging over x\ we get (e/L)v which is essentially 
the velocity operator. 

In the absence of driving the "pure" stationary states of the system are the 
microcanonical states. We use classical language but also have in mind a semiclassical 
picture. Each microcanonical state occupy a shell whose phase space volume is {2-Kh) d . 
The density of states is 

The zero order stationary state is characterized by an occupation function f(E). Later 
we shall take it to be the Fermi function. Thus 

(12) 

where E = H(r,p) is the energy. The distribution of the particles in energy is 

p(E) = g(E)f(E) (13) 
and the total number of particles is 

N = j j p(r,p)^ 1 = J~j(E)dE = J°°J(E)g{E)dE (14) 
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3. Network modeling 

A network is defined as a set of ID wires that are connected in vertices. The network 
Hamiltonian is ill defined in the classical limit because the the scattering in each vertex 
is described by a scattering matrix. In particular for the model of Fig. Id the scattering 
is described by a scattering matrix S a b, and we define the corresponding transition 
matrix as g a b — |Safc| 2 - Thus the classical description of the system is stochastic rather 
than deterministic. 

Still we can regard networks as an effective way to describe the chaotic dynamics 
|19j . The reason is that upon coarse graining a chaotic system looks like a stochastic 
model. Specifically in the case of the system of Figlb, quantum mechanics introduces 
"coarse graining" in a most natural way. Each mode in the scattering problem can be 
regarded as a ID wire with the dispersion relation 

V ^if^-Wf (15) 

where E is the energy of the particle, and L± is the width of the waveguide. The open 
modes are are those for which v n is a real number. We shall denote the number of 
open modes by M. hence the number of open channels in the open geometry is 1M.. 

The density of states of the system can be written as a sum over single-mode 
expressions: 

M 

g(E) = g dot {E) + 2 £ (16) 
^ (27rft) v n (E) 

where the factor of two takes into account both clockwise and anticlockwise motion. 
A stationary state of the system is described by the distribution functions p^(E) of 
the clockwise moving particles and p^~(E) of the counter-clockwise moving particles. 
The index n distinguishes different modes. The normalization is such that 

N= f Pdot (E)dE + LJ2 I (P7(E) + p7(E))dE (17) 

The density of particles per unit length in a given mode is implied by Eq. l|13|l : 

dE 

P ^{E)dE = P ^(E)dE = f(E) (18) 

{2Trh)v n (E) 

Note that for a microcanonical distribution dE can be regarded as a fixed parameter 
that defines an energy window or a width of an energy shell. 

The scattering is described by a 2Ai x 2A4 transition matrix that has a block 
structure: 

gab= { 9 g T f 11 ) (19) 

It consists of the reflection matrix g^ m and the transmission matrix g^ m . Note that the 
channel index a contains both mode specification and left /right lead specification. We 
assume time reversal invariance, so as to have a symmetric matrix. For clarification we 
note that if N particles incident in channel b, then g a bN particles emerge in channel a. 
This means that g a b is the ratio between ingoing and outgoing fluxes. In the ergodic 
state Ea. 1)18(1 implies that p a cx l/v a . Therefore we have detailed balance: 

g a b PbVb = gba p a v a [no summation] (20) 
Namely, for a stationary state the transitions from a to & are exactly balanced by the 
transitions from b to a. 
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4. Definition of the conductance 

If the EMF is not too large we expect to have Ohmic behavior (T) oc EMF. Hence we 
define the conductance via the relation 

(I) = -G$ (21) 

This is in fact a special case of a more general concept of conductance. In the theory 
of "quantum pumping" we may vary in time some other parameter X, and get current 
(X) oc X. Hence we can define (generalized) conductance via 

(I) = -GX (22) 

Much of the formalism that we are going to use can be extended to handle this more 
general case. However, in what follows we restrict ourselves to the case of an EMF 
driven system. 

Still one can wonder whether it is important to specify how the EMF is introduced. 
Does it matter how a(x) look like? The answer is that within linear response theory 
all the choices lead to the same result. This is not merely a gauge issue because for 
different a{x) the electric field along the ring does not look the same. Namely, if the 
gauge of the vector potential is changed 

a(x) := a(x) + (d/dx)v{x) (23) 

then the new scalar potential is 

V(x) := V(x) - e<S>v(x) (24) 

Thus a different choice for xq can be regarded as associated with adding a rectangular 
barrier of height e^. Within linear response it is assumed that e<I> is too small to 
make any difference. If this assumption is not applicable, it is no longer the "linear 
response regime", and the specification of a(x) becomes significant. 

5. The long time scenario 

Some people find it inappropriate to define conductance for a closed system because 
the problem does not possess a stationary solution |22| . Namely, it is clear that without 
a contact to a thermal bath the driven system is gradually heated up. However, we 
find this objection of no relevance. The practical point of view of an electrical engineer 
is demonstrated in Fig.lc. It is clear that at any moment the engineer is inclined to 
characterize the ring by its conductance. This is true irrespective of whether there is 
a contact with a thermal bath or not. In the absence of such contact it is evident that 
the system is heated up and therefore the conductance becomes time-dependent. 

It is true that the overall scenario is (formally) beyond linear response, but it is 
also true that at a given instant of time it is feasible to have a valid linear response 
description. The validity condition is having an ergodization time which is much 
smaller than the time that it takes to have a significant change in the (evolving) 
energy distribution of the system. This reasoning leads to "slowness conditions" that 
are further discussed in Ref . [T7| . 

We note that there is a strict analogy here with the strategy of analyzing quantum 
pumping. Also there the conductance G(X) is calculated, using either BPT or the 
Kubo formula, for various points X — (X±, X2) in parameter space. Later the pumped 
charge is expressed as a line integral over the conductance: 

Q = (f(I)dt = <b[-G(X(t)) ■ X]dt = - fG-dX (25) 
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Consequently the pumped charge is not proportional to the amplitude of the driving 
cycle. Thus we can get (globally) "non linear response" out of (momentary) linear 
response analysis. See jS] for more details. 



6. Digression: The role of the environment 

On the mathematical side we have already defined precisely the assumptions in the 
basis of the derivation that we are going to present. The purpose of the present section 
is to further clarify the physical circumstances that justify these assumptions. Also 
we would like to put the present work in context of past literature. This section can 
be skipped in first reading. 



6.1. Quantum to classical correspondence 

Our first step is to define a feasibility condition for the validity of a semiclassical 
treatment. Assuming that the total transmission of the device is gr ~ 1 the time for 
the randomization of the velocity is 

r cl « (r^—) x — (26) 

where L is the length of the ring, and Vf is the Fermi velocity. On the other hand the 
time that it takes to resolve the quantized energy levels of the ring is 

iH 0i3 e„ b0IE « M x — (27) 

VF 

Hence the quantum-to-classical correspondence condition is 

M > — — (28) 

i-gr 

which is always satisfied in the classical limit. Note that the limit M. — > oo is analogous 
to h->Q. 



6.2. Single mode rings 

In the case of a single mode ring (a one dimensional ring with a delta scatterer) 
a classical treatment of the dynamics does not make any sense in view of the 
correspondence condition of the previous subsection. The analysis should be 
purely quantum mechanical and issues such as Landau-Zener transitions |12j . Debye 
relaxation mechanism [Jj and Dynamical Localization E] should be taken into 
account. Recently we have introduced a new ingredient into the analysis that sheds a 
new light on the whole issue 14 . 



6.3. Multi mode diffusive rings 

In the case of multi mode diffusive ring the leading order result for the conductance 
is as expected just the classical Drude expression. The typical calculation [H] assumes 
that the levels are "broadened" due to the interaction with the environment. In the 
major case of interest the level broadening T is assumed to be larger than the mean 
level spacing A but much smaller than any semiclassical energy scale. Hence it barely 
affects the Drude result. Still it determines the quantum weak localization correction, 
which turns out to be of order A/T. The weak localization correction depends on the 
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levels statistics and therefore on the magnetic flux. Hence it can be detected in an 
actual experiment 

6.4- Multi mode ballistic rings 

In the present work we consider a multi mode ballistic ring rather than a diffusive 
ring. This means that the time to randomize the velocity can be much larger than the 
time L/vf to make one round along the ring. This implies that in our configuration 
the conductance (in natural units) can be larger than the number of open modes. 

We assume that the environment induces level broadening T which is larger than 
Hvf/L. Hence we derive the leading (classical) term and do not take into account the 
implications of quantum interference. To put this assumption in a larger perspective 
we make the following classification: 

f . Isolated system (no environment). 

2. The bath induces only decoherence effect. 

3. The bath induces velocity randomization. 

4. Bath limited dynamics. 

The first case of fully coherent dynamics will be analyzed in ^5] . We consider in this 
work only the second case. Still we would like to further explain why the second case 
is physically typical, and to make some comments on the other two cases. 

6.5. The decoherence mechanism 

It is typical to assume that the fluctuations of the environment are of large spatial 
correlation length compared with L±. An extreme case is the Caldeira-Leggett 
modeling which assumes an infinite correlation length. The matrix elements of 
the position variable scale like L± for inter-mode transitions and like L for intra- 
mode transitions. Hence inter-mode transitions are rare compared with intra- mode 
transitions, the ratio being (L±/L) 2 . Therefore it is realistic to consider circumstances 
such that velocity randomization due to the environmental "noise" can be neglected, 
while intra-mode transitions cannot be neglected. The latter lead to decoherence. The 
simplest estimate for the decoherence rate is T = r/ksTL 2 /h, where r\ characterizes 
the coupling to the environmental modes, and T is the temperature. Thus we see that 
case 2 is physically typical. 

6.6. Velocity randomization 

If the fluctuations are strong enough inter-mode transitions cannot be ignored. This 
would lead to randomization of the velocity in the wire region. The scenario of having 
a bath that just randomizes the velocity, but does not affect the transmission of the 
ring is apparently not of much physical interest. Still such effect can be realized 
artificially, and it is of pedagogical importance, as discussed in section 12. 

6. 7. Bath limited dynamics 

If the interaction with the environment determines the transmission of the ring, we 
get to case 4. The most obvious example is the the analysis of conductance in room 
temperatures. The scattering and hence the diffusion of the particles is dominated by 
bath induced inelastic scattering by phonons. 
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The main point regarding case 4 is that the bath cannot be eliminated from the 
model analysis. Another example for such circumstances is provided by the Debye 
dissipation mechanism [7|- The latter assumes that the inelastic relaxation time is 
much shorter compared with the adiabatic variation of the energy levels, leading to 
time lag between the driving cycle and the adjustment of the occupation probabilities. 



7. Linear response analysis 

The simplest route to linear response theory |121 1161 1171 ITS) takes the relation 
dH/dt = dH/dt as a starting point. It follows that the change in the energy of a 
particle is given with no approximation by the formula 

H(t) - H(0) = -$ / T{t')dt' (29) 



By squaring and averaging over initial conditions we get that the second moment as 
a double time integral over (J r (t')J r (t")). Within linear response this correlation is 
approximated by the stationary correlation function 

C{t' -t") = (T{t')T{t")) E (30) 

where the average on the right hand side is taken with a zero order microcanonical 
solution. Thus one concludes that there is a diffusion in energy, with the coefficient 

1 f°° 

D E = $> 2 x- C(r)dT (31) 

* J -oo 

Next one wants to see what happens in the more general case of an arbitrary f(E). 
On long times it is argued that the probability distribution p(E) satisfies the following 
diffusion equation: 

l = 4U4f^) (32) 



dt dE V dE \g{E)' 

The energy of the system is (H) = J Ep(E)dE. It follows that the rate of energy 
absorption is 

d_ 

dt v *' " Jo ~~ ~ a 8E \g{E)J 

For zero temperature Fermi occupation we get 
d 



-(H) -J™dEg(E) D E A (ztS) ' ; «) 



dt 



(H) = g{E)D E = (34) 



E—Ef 



This is the mesoscopic version of Joule law. The expression for the conductance is 
1 
2 



G = g(E F ) x 1 / C(r)dT (35) 



8. The Kubo formula 



The above is apparently the simplest and most illuminating derivation of the 
Fluctuation-Dissipation version of the Kubo formula. A more complicated treatment 
[T%1 0] 03 allows to write a generalized version that holds also for "quantum pumping" 
applications. Namely, 

G = g(E F )x / C{ T )dr (36) 
Jo 
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with 

G(r) = {T{t)T{Q)) e (37) 

where T is the generalized force which is associated with the driving. In the present 
application T is just the current operator, hence C(r) is symmetric, and therefore the 
generalized version Ea. l|36[) is equivalent to Ea. (|35[) . 

It is also important to make a connection with the more traditional treatment of 
conductance in case of disordered metals. If we set a(x) = 1/L for the vector potential 
we get (e/L)v as the current operator. Hence we get from Ea. (|35|) 

G = g(E F )x-(-) / (v(r)v(0))dT (38) 



2 \L 

The conventional derivation of the Drude formula is based on the assumption of 
exponential decay of the velocity-velocity correlation function. The latter formula 
is (formally) valid in any case, but in case of a disordered sample it has special appeal 
because it implies the Einstein relation between the conductance and the spatial 
diffusion. Namely, if we have a diffusive ring then it is natural to write G = G/L 
and to define g(E F ) = g{E F )/L. Then we can rewrite Eq.|j3HJ) as 

G = e 2 g(E F )D spaco (39) 

The more conventional derivations of this expression is based on the phenomenological 
relation J = —DV (density) — GV (potential) and the argument that J = at 
equilibrium. 

9. Implications of the Kubo formula 

From the the Kubo formula Ea. H36|) it is not obvious that the result for G is 
independent of where we measure X. In general it can be proved |20| that for a different 
choice of T the corresponding Kubo conductance may differ at most by eg(E F )(!F). 
But if T is a current operator then (T) = and therefore the conductance be becomes 
independent of X\ . 

Also it seems that the Kubo conductance is proportional to the density of 
states. Therefore, if we had doubled the volume of the cavity, would we get larger 
conductance? Furthermore, does the result for G depends merely on the transition 
matrix, and not (say) on the dwell time inside the scattering region? 

To answer these questions, and to establish the x\ independence of G, let us write 
the Kubo formula in a more illuminating way. By definition we have 

/oo 
(T(r)^(0))dT = YtPr^rQr 
-oo r 

= Wjl 5{E ~ n (r, P ))^F(r,p)Q(r,p) (40) 

where r is an index that labels phase space cells (different initial conditions), and p r 
corresponds to a microcanonical distribution. We have introduced the notation 

/oo 
I(T;r,p)dT (41) 
-oo 

Namely Q(r,p) is the total charge which is obtained by integrating the current which 
is induced by a particle that goes through the point (r,p) at t = 0. It is in fact 
(for e = 1) the winding number of the associated trajectory, and therefore it gives a 
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result which is independent of the chosen section. Note however that Q(r,p) obtains 
a meaningful value only upon course graining, else it is erratic. Now we can write the 
Kubo formula as 

G=IJS(E- H(r,p))^F(r,p)Q(r,p) (42) 

This expression has several advantages. One advantage we have mentioned: the 
result is manifestly independent of the definition of the current operator. The second 
advantage is that it shows that the global density of states is in fact not important. 
We can double the volume of the cavity, still we would get the same result provided 
that the scattering probabilities are not affected. In particular we see that time delays 
are not important. 



10. Conductance of a single mode ring (part 1) 

In this section we show how the Kubo formula for a closed ring leads to a 
Landauer-alike formula for the conductance provided the effect of the environment 
is to completely randomize the velocity within the wire region without affecting its 
transmission. To simplify the presentation we consider the single mode case. 

For the Kubo formula we have to evaluate the correlation function C(r) of Eci. (|37|l . 
and to calculate the integral in Ea. (|36|) . As explained the result of the calculation 
should be independent of the how we define the current operators. The simplest choice 
is to define T as the current through a section x = xq on the left of the scattering 
region, while X is the current through a section x = x\ on the right of the scattering 
region. C(r) comes out as a sum of delta functions. The shortest time correlation is 
associated with the time t± to cross the scattering region. For example, if there is no 
time delay then t\ = (x% — xq)/ve- We can regularize T as a rectangular of width 
e. The probability to have there a particle moving in the right direction, such that 
JF(0) = evE/s is (e/L)/2. The current that we get in the other side of the barrier is 
2 = eS(r — ti). Assuming that this is the only correlation, and taking into account 
the time reversed correlation for r < 0, we get 

Note that if we had chosen x\ — x$ we would get three delta functions: a self 
correlation delta function 5(t) and reflection peaks. Namely, 



C(r)=e^ 



5( T )-i(l- 5T )^5(r T To) (44) 
z ± 

where To is the scattering time. Obviously the integral over the new C (t) is the same 
as the integral over the former one. Irrespective of our choices we get from Kubo 

G = &H 9T (45) 
which looks like the (single mode) Landauer formula. 
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11. Conductance of a single mode ring (part 2) 



If the velocity is not randomized within the wire, then there are other correlations 
that involve the time L/ve to encircle the ring. For the following calculation it is 
simplest to assume that < xq < X\. What we have is to calculate the integral 



{l{ T )T{Q))dT = Y^ Vr T r Q r 



(46) 



Note that we find it convenient here to set r = as the lower bound of the integral. 
Recall that r is an index that labels phase space cells (different initial conditions), and 
p r corresponds to a microcanonical distribution. Note that T is non-zero only if r is 
located at x = xq. The total charge which is transported through the section x = x\ 
is defined here as 



(X(t)Mt 



(47) 



where the current is evaluated under the assumption that the particle is launched at 
point r. There are two relevant possibilities: Either the particle is launched at x = xq 
in the clockwise direction, or it is launched at x — x$ in the anti-clockwise direction. 
Observe that the (net) charge that goes through the section x — x\ after a round trip 
is suppressed by a factor (2<7t — 1) due to the scattering (we sum the clockwise and 
the anticlockwise contributions). Thus we get that the total charge that goes through 
the section is 

= e [1 + (2g T - 1) + (2g T - l) 2 + ...] = 1 



2(1 -g T ) 



(48) 



for a particle that is launched clockwise, and 
1 1 



.2(1 -gr) 

for a particle that is launched anti-clockwise Thus we get 
\l(T)F(0))dT = ±.(+ev E )Q- + -L(_ eVB )Q- = 



L 



9t 



9T 



leading to 



G 



e 

2irh 



9t 



9T 



(49) 



(50) 



(51) 



12. Relation to the Landauer formula 



As explained in section 10 we get from the Kubo formula a Landauer look-alike formula 
if we assume that the environment induces velocity randomization in the wire region 
without affecting its transmission. In fact we can get to the same conclusion by 
modeling the "loss of memory" in the wire region as a scatterer with transmission 
g wiro = 1/2. It is well know that the G of Ea. (|51|) obeys Ohm law for addition of 
resistors in series. Hence 

-l 



G 



e 

2irh 



9t 



1- 



fjT 



e 

2~Trh 



9o 



1-9, 



2irh 



9o 



(52) 



We would like to emphasize that the purpose of this section is purely pedagogical. 
As stated in section 6 an environment that just randomize the velocity without 
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affecting the transmission is apparently of no physical interest. Still if one insists 
it can be constructed artificially. Simply cut the wire and connect the two ends to a 
chaotic cavity. A particle that moves in the wire gets into the cavity and after a time 
delay gets out either clockwise or anti-clockwise with equal probabilities. Hence in 
such arrangement g vire = 1/2. 

The pedagogical importance of the above discussion is in making a bridge between 
the reservoir philosophy of the Landauer construction and the Kubo formalism of 
closed systems. The memory loss device that we have described above provides the 
same "service" as the reservoirs in the Landauer picture. 



13. Conductance of a multi-mode ring 



Let us assume that the EMF — $ is concentrated across the scattering region (x = 0). 
When a particle goes through x = it gains momentum (p t— > p — e$/w), where 
v = |p|/m. Hence the change is energy is E £^ e$ for right and left movers 
respectively. The state of the system is described by the distribution functions p^ (E) 
and p^(E) of Ea. 1)17(1 . The index n distinguishes different modes. It is implicit from 
now on that we look for an ergodic-like solution, such that the density of the particles 
along the ring is uniform. The balance equations are: 



dt 
dt 



= - [pn v n] 



[Pn v n 



y 

9nmPm Tfi 



y 



^ t 9nmPm v m 



y 

9nmPm 



(53) 



(54) 



It can be verified that the zero order (<j> = 0) stationary solution of this equation is 
given by Ea. (|18fl . where f(E) is an arbitrary function. We are looking for a first order 
stationary-like solution. The linearized equation for the clockwise moving particles is 



E^ ~ 9 T )n m Sp^v n 



y9nm S PZVn 



2ir dE ? ^' 



(55) 



A similar equation exist for the counter-clockwise particles, 
corresponding equations we get 



/ j ilnra 
m 

Subtracting the 



E 



[(1 - g T ) + g R }„ m [p^v - p^v] m = 2^$M^ 5 „ 

2-7T dE 



with the solution 

[p^v - p^v] 

The current is 

* = E 



= $— 

2tt 



2ir dE 



(l-3 T )+.9 i 



(56) 



(57) 



dE (Pn - Pn)ev n 
2 



E 

nm 



i 



g T + g R 



df(E) 
dE 



dE 



(58) 



With the assumption of Fermi occupation we get Eq.©. Note that upon summation 
the order of matrix multiplication is not important because g nm is symmetric. 
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14. The wire with cavity model system 

We consider a ring (Fig. lb) which is formed by folding a rectangular waveguide (i.e. 
imposing periodic boundary conditions). A chaotic cavity is attached to the waveguide 
at one "point" . A particle has some probability to enter the cavity, where memory is 
"lost" , and then it gets out again with equal probability to either side. A particle that 
travels in mode n of the waveguide has a transverse momentum ±(irh/L±)n where 
L± is the width of the waveguide. The distance between subsequent hits of the same 
wall is 

y/2mE- ((Trh/L ± )n) 2 ^ r VM 2 



step 



-2Li = 



Tl 



2Lj 



(59) 



(irh/L±)n 

The number of open modes M. is implicitly defined via the latter equality. The 
probability to get into the cavity via an opening of size L op is: 



Pn 



step 



minimum 



1 



(60) 



where a = L op /(2L±). The crossover from p n < 1 to p n = 1 happens at 



n c = 



(61) 



We are going to treat M. as a free parameter. Hence we have two parameters that 
characterize the scattering: the classical (geometrical) parameter a, and the quantum- 
mechanical parameter A4. Note that the classical limit is A4 — > oo. 

Let q n be the probability to get out of the box to mode n, either to the right 
going channel or to the left going channel. It follows that g^ m = (l/2)q n p m . From 
9nm = 9mn we conclude that q n /p n = c is the same for all channels. Taking into 
account that ^2 n q n = 1 we get c = l/(J2 n p n ), an d hence 
1 



Q R 

cJnm 



T 
Qnm 



q-nPn 



TPnPr, 



+ (1 -Pm)S 



= -jPnPm + (1 - Pm)S r , 



(62) 
(63) 



Thus, given the input parameters a and A4, we can calculate g nm . It is useful to 
define the total probability of transmission for a particle that comes in channel n as: 

1 



9n 



= i 



-,Pn 



(64) 



For sake of later estimates we note that for M. 3> 1 , sums over n can be approximated 
by an integral over x = n/ Ad. Using the obvious notation x c — n c /M we get 



w af\{l/xf-l]-^dx+[l-. 



1 

a 2 



l + -a- 



(1 



,2^1/2 



+ [1 - X C ] 



(l + a 2 )-V2 



1 + -0- 



(65) 
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and 



1 

M 



l/M 
111 

111 



[(1/x) 2 ~l] 1/2 dx 
2x c 



^2)1/2 



(1 



a 2)l/2 



M 



M 



(l-(l-^)Va) 



[ln(2M) -11- 



0(a) 



(66) 



We now turn to the calculation of the conductance. First of all, let us calculate 



the Landauer conductance. Thanks to the simple structure of the g\ 
calculation is quite easy 



e 

2^h 



'7im 



— Y 



9n 



matrix, the 



(67) 



Each channel has total transmission in the range (1/2) < g n < 1 and therefore 
the conductance (in normalized units) roughly equals to the number of open modes. 
Substitution of leads to: 



— T(i 



(p»/ 2 » = hi [1 -i a+<D{a4)]M 



(68) 



For the multimode conductance of Eq.@ the calculation is more complicated. 
At first sight it seems that the calculation should be done numerically as in Fig. 2. 
The numerical calculation in Fig. 2 (circles) is done in a way which is inspired by 
a similar type of calculation within the Landauer formalism. We define h nm = 
([2g T /(l — g T + /J]™) 1 / 2 and write the sum in Eq.© as trace [h^h] . Then we 
make singular value decomposition of h and sum over the squares of its eigenvalues. 

The other way to calculate the multimode conductance starts with an attempt 
to make a zero order evaluation of the sum. This means setting c = in Ea. (|63|) . 
The resulting estimate gives a rough approximation as seen from Fig. 2 (crosses). The 
main source of error are evidently the low modes. Surprisingly it turns out that the 
calculation can be carried out to infinite order in c, thanks to miraculous cancellations. 
Using the expansion 



1 



A-cB 



1 11 

A A A 



o 1 1 1 

c 2 -B-B- 
A A A 



(69) 




200 



Figure 2. The left panel displays p n for a system with a = 0.7 and M = 200 open modes. 
The crosses in the right panel are for the cumulative sum over (1 — pn)/pn- The circles are 
calculated from the exact formula. Namely, the matrix h is diagonalized, and the cumulative 
trace over the square of its elements is displayed. 
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with A nm =p n S nm and B nm = (\/2)p n p m , we get the result 

= Snm + (c/2) PnPm h 

m Pn Pn Pm 



(c/2) 2 V] —PnPk—PkPni — 



Pn 



Pk 



Pn 



= S nr . 

Pn 



(70) 



where in the last step we have made a geometric summation over all orders. Now we 
can calculate the conductance 



G 



-T. 



9 



1 



9t 



e 

2TTh 




cM - M 



(71) 



Recall that c = l/^2 n Pn- Hence this expression requires merely the evaluation of 
the sums ~}2, n Pn and If we have p n » 1 for all modes, then we get simply 

G = (e 2 /(27rfr)).M which reflects that number of open modes. But the interesting case 
is when a is small: 



G 



e 

2^h 



(1 + hx(2M)) - 1 + 0(a) 



M 



(72) 



Unlike the case of the Landauer conductance, the result does not reflect the number 
of open modes. The contribution of the low modes is singular in the limit of small a. 
Furthermore, the conductivity (conductance per channel) diverges logarithmically in 
the classical limit. 



15. Concluding remarks 



Much of the derivations in this paper can be generalized to analyze "quantum 
pumping". In Ref.|21| we have considered a single mode device where the current 
is induced by translating a scatterer. Namely, X — —GX and hence the transported 
charge is dQ — —GdX where dX is the displacement of the scatterer. From the BPT 
formula one obtains 



G 



-(1 - g ) x -k F 



(73) 



where go is the transmission of the scatterer and k F is the Fermi momentum. If we 
close the system into a ring, and use the same assumptions as in this paper we get 

1 - go' 



G = — 



9o 



9t 



1 



gr 



(74) 



where gr is the overall transmission of the device. Multi-mode generalizations of these 
results can be obtained by employing either the Kubo or the Master equation approach 
as in the present paper (not published). 

We regard the "Ohmic" problem that has been discussed in the present paper, 
and the above mentioned "pumping" problem of Ref.|21|. as the prototype models 
for the application of linear response theory: The "Ohmic" problem has to do with 
the dissipative part of the response, while the "pumping" problem has to do with the 
geometric (non-dissipative) part of the conductance matrix. In both cases we can use 
the Kubo formalism as a starting point, and in both cases we can regard the scattering 
approach (Landauer-BPT) as a special limiting case. However, one should be aware of 
the subtle differences between the two problems. The main point to remember is that 
the adiabatic limit of "quantum pumping" is the non-vanishing "adiabatic transport" 
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formalism, while the adiabatic limit of the "Ohmic" problem gives zero conductance. 
For further details see 

We have assumed that the coherence time is short compared with the time that it 
takes to encircle the ring. This makes the calculation insensitive to the energy E. Once 
we consider a strictly isolated system, we get an E dependence that has to do with 
semiclassical energy scales such as (1%/L)vf- Such energy scales are larger compared 
with the mean level spacing (oc h d ), and may invalidate the Kubo formalism. An 
extreme example for the implication of having a non-universal energy scale is analyzed 
in in Ref.JI]. The study of the general multi-mode case introduces further conceptual 
as well as technical complications |15j . 
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